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Abstract. The cosmological constant A modifies certain properties of large 
astrophysical rotating configurations with ellipsoidal geometries, provided the objects 
are not too compact. Assuming an equilibrium configuration and so using the tensor 
virial equation with A we explore several equilibrium properties of homogeneous 
rotating ellipsoids. One shows that the bifurcation point, which in the oblate case 
distinguishes the Maclaurin ellipsoid from the Jacobi ellipsoid, is sensitive to the 
cosmological constant. Adding to that, the cosmological constant allows triaxial 
configurations of equilibrium rotating the minor axis as solutions of the virial equations. 
The significance of the result lies in the fact that minor axis rotation is indeed found 
in nature. Being impossible for the oblate case, it is permissible for prolate geometries, 
with A zero and positive. For the triaxial case, however, an equilibrium solution is 
found only for non-zero positive A. Finally, we solve the tensor virial equation for the 
angular velocity and display special effects of the cosmological constant there. 
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1. Introduction 

In the last decades the universe has been mapped by different methods sufficiently 
accurate to allow us global and detailed conclusions about its properties and behavior. 
Key observations are the luminosity-redshift relationship of supernovae of type la (SNIa) 
PP, the matter power spectrum of large scale structure as inferred from galaxy redshift 
surveys [2|, and the anisotropies in the Cosmic Microwave Background Radiation 

These observations allows us to conclude that the Universe is spatially flat, and 
surprisingly its expansion is accelerating (instead of the long thought slowing down). 
Many different models have been proposed to explain this acceleration: For instance, 
large scale modifications of General Relativity 0]. The most popular explanation is, 
however, the inclusion of a dark energy component which corresponds to about 70% of 
the total energy budget of the universe and has a quite negative pressure. Candidates 
for the dark energy range from scalar fields [3] (which are widely used in cosmology [S| ) , 
up to dark energy/dark matter unifying candidates, as for instance the Chaplygin Gas 
[7j. The simplest model, however, would still be a positive cosmological constant A. 

Understanding the nature of dark energy is one of the major objectives in cosmology 
nowadays. One way to undisclose its properties is to investigate its astrophysical and 
cosmological imprints. Dark energy affects the formation of galaxies and clusters [HUH], 
as well as dark matter halo properties, such as the virial mass and the virial radius 
In ^T] it was also shown that A affects the motion of test bodies in the Schwarzschild-de 
Sitter metric over large distances, a case of importance for large scale objects. Although 
observables within the solar system seem not to be sensitive to the global expansion of 
the universe that is not the case with large small-density astrophysical objects. In 
fact, investigations on the effects of dark energy in the local Hubble flow show that those 
effects may not be negligible [T^j. In case of a cosmological constant, it has been shown 
that not only the formation of astrophysical bodies gets affected by A, but also its static 
properties like the internal mean velocity and the angular velocity ^3]. In fact, it was 
shown in ^H] that A may produce lower values of galaxy velocity dispersions. Taken all 
these findings together one can state that dark energy has some effects at astrophysical 
scales. It is then important to try to understand what further possible effects could dark 
energy imprint on some of the properties of the astrophysical bodies. These signatures 
can be used to probe and test dark energy candidates. 

The effect of A is bigger the more we deviate the object's geometry from the 
spherical case fTJ • This is important since by far not all structures found in 
the universe have spherical symmetry. Rotating configurations acquire non spherical 
geometries which arise as solutions of the tensor virial theorem (of second order). 
With zero cosmological constant some of these figures of equilibrium correspond to 
Jacobi ellipsoid (triaxial, rotating along the minor axis) and Maclaurin ellipsoid (oblate, 
rotating along the minor axis). Also one finds the Dedekin ellipsoid, (triaxial and flat 
because internal streaming), the Riemann ellipsoid (triaxial with rotation and internal 
motion), among other possible solutions. Interestingly, not only galaxies have different 
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morphologies [TH], but also clusters of galaxies. Of interest for us is especially the fact 
that clusters can have a prolate shape (hence we expect that they can also appear 
in triaxial configurations) and they can rotate j!8j . 

In this article we will investigate some basic properties of ellipsoids in the presence of 
a cosmological constant. In particular, we will answer questions related to the geometry 
of a rotating configuration. For instance, given the axis of rotation, which could be 
any one of the three axis, what is the allowed range of the axis which emerges as a 
solution? A legitimate question is also what are the changes regarding the bifurcation 
points? In other words, at which point is a triaxial solution allowed. Finally, we can ask 
what is the possible angular momentum? Although not frequently encountered, there 
are astrophysical objects rotating along the major axis (this is the same to say that 
the objects rotates the minor axis) (201 EI]- It is interesting to make some theoretical 
attempts to understand such exotic configurations. We will show that the triaxial case 
[23] of minor axis rotation is possible with positive A and the allowed region associated 
to values characterizing the geometrical configuration vanishes with vanishing A. For 
oblate and prolate case we calculate the bifurcation points emphasizing the modifications 
and new effects induced by A. Clearly, the ellipsoidal geometry which we will study in 
the paper is well suited to model many geometries such as the disk. Therefore, the 
results of this article should in principle be applicable also to superclusters which are 
usually found to be very flat objects [TH] . 

Hereby we will concentrate on homogeneous rotating ellipsoids, whose equilibrium 
configuration will be explored by the tensor virial equation with A [2E1 El] describing 
the equilibrium properties of possible non-spherical configurations. The tensor virial 
theorem results into two equations: an equation for the angular velocity and one 
restricting possible geometries. It is therefore a perfect tool to study the subject. It is 
known that for spherical systems it is sufficient to work with the scalar virial theorem. 
Our results generalize the ones obtained in [2BJ- This is especially true for the bifurcation 
points of rotating ellipsoids as solutions of tensor virial equations for the oblate geometry. 
Indeed, we follow closely the exposition and methods given in |2H] which applies for the 
case A = 0. However, the study of the prolate/triaxial case with minor axis rotation is 
new. It is not surprising after all that this has not been examined before as this case is 
only possible with A ^ 0. 

This paper is organized as follows: in section 2, we review the basic expressions 
for describing equilibrium in the Newtonian limit with a cosmological constant. Since 
we retain in this limit A, we are essentially dealing with a non-relativistic de Sitter 
spacetime (also called Newton- Hooke). Section 3 is devoted to a brief discussion about 
the applicability of our approach to galaxy clusters. In section 4 we explore the allowed 
configurations for a rotating ellipsoid whose rotation lies along the minor axis (major 
axis rotation). The next section explores the strange scenario, which corresponds to the 
case of minor axis rotation, i.e. when the angular velocity is along the major axis. We 
end with conclusions. 
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2. The Newtonian limit with cosmological constant 

To explore the consequences of a positive cosmological constant, we make use of the 
Newtonian limit of Einstein field equations. In this limit, the Poisson equation is written 

as 

V 2 $ = 4nG N (p + 3^) - A. (1) 

where P is the pressure, which arises as a general relativity correction. The solution of 
((TJ at the zeroth order ofv/c (from now on we set c = 1) takes the form 

Q{r) = -G N f ^-dV-iA|r| 2 , (2) 
Jv |r — r'| o 

where we have neglected secondary effects of the cosmological constant in the set of the 
Dirichlet boundary conditions (2H1- The potential is then a contribution from a pure 
gravitational part $ grav and a contribution from the expansion proportional to A. 

Together with equation (JTJ), a self gravitating configuration is ruled by an equation 
of state and the hydrostatic equilibrium equation (Euler's equation), which is written 
as 

dll ' 

p —l +pd .<S> + d . V .. = : (3) 

where Vij is the pressure tensor (which for simplicity we take as isotropic Vij = SijP) 
and Ui = (vi) is the mean values of the velocity. The second order tensor virial equation 
can be derived by taking spatial moments in Euler's equation. The result is an integral 
equation 



ld 2 X 



2T ik - | Wife | + -vrGivp vac Xifc + U ik , (4) 



2 dt 2 1 1 3 



with 



n ifc = / V ik d\ T ik = - [ p(vi)(v k )d 3 r, (5) 
Jv 2 Jv 

Wk = -G N f p(v) n d k ^dh. (6) 



which corresponds to the pressure integral, the kinetic energy tensor and the 
gravitational potential energy tensor, respectively. The moment of inertia tensor is 
defined as 



J prir k d 3 r. (7) 

In statistical mechanics the virial theorem follows from taking time averages over 
a long period of time (denoted below by (...)) such that the configuration reaches 
equilibrium. This implies that in equilibrium, the left hand side of (J3J) vanishes which 
gives us the tensor virial theorem as 

2(T ik ) T - (|y\4|) T + ^G NPvac (l ik ) T + (U lk ) T = 0. (8) 

In astrophysics this would imply that we have to take an average over ensembles of 
observed objects of a specific type. This, in principle, is possible noting that the fact 
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that the external contribution, i.e, the vacuum energy density p vac = A/87T is constant, 
and therefore the time average is only taken over proper parameters of the configuration. 
However, it is also customary and more convenient to assume that for most astrophysical 
objects the second derivative of the inertial tensor vanished at a sufficiently large time. 
The averaging is not necessary then. For generalizations (dark energy , scalar fields or 
Chaplygin gas), where one has also to take into account the evolution of the background, 
this is a more straightforward approach. 

In the next sections, we drop the (...) brackets and also use the convention 
G N = c = 1. 

3. The virial theorem and large astrophysical structures 

It is of some importance to this work to discuss briefly the question whether the virial 
theorem or hydrostatic equilibrium apply to large astrophysical structures. This is of 
importance for the present study because certainly the effect of A will be seen only for 
relatively low density objects, i.e. clusters and superclusters. The other issues which 
are of phenomenological interest for us is the shape of these objects and their possible 
rotation. 

Whether one can apply the virial theorem (or hydrostatic equilibrium) and its 
consequences to groups or clusters of galaxies has been a matter of debate over years 
[3H1IHUIII2]- In many cases (and these are the cases to which we can apply our subsequent 
results), this dispute has been decided in favour of the virial theorem [3TJIIIII]- Indeed, 
over decades the virial technique is the most widely used one |SH ES| and the agreement 
with other methods is satisfactory [^H] • The assumption of hydrostatic equilibrium in 
clusters is common and, as it seems, works well jHZj • There are exceptions due to cooling 
flows or misidentification of the cluster shape (i.e. superposition of two apparent 
individual objects) [39;. Last but not least, the virial theorem with cosmological constant 
has been successfully applied to galaxy clusters at early epochs |HJ- 

The Morphology of clusters is as rich as the corresponding one for galaxies |16j . 
Almost all shapes can be fitted into ellipsoids [13 CHI HOI HI] • Interestingly, virial results 
are used to take into account finite size effects of galaxy clusters in lensing |42j . These 
effects are parameterized by ellipticity, a concept related to eccentricity of ellipsoids. 
The second important feature of clusters, mentioned already in the introduction, is that 
they can exhibit rotation JH] exactly as galaxies do. The applicability of virial theorem 
to clusters taken together with the facts of their ellipsoidal nature and the established 
rotation means that we can apply our results below to galaxy clusters. It remains to 
address the question whether the application is valid for the luminous matter or the 
'whole' object including dark matter. We can rely here on studies and simulations 
which all, and this is the crucial observation here, indicate that dark matter halos 
around clusters are of ellipsoidal nature [HH EI] • However, we could apply the virial 
theorem to the luminous structure alone including for the total mass also the dark 
matter component. This seems possible as the form of the luminous matter does not 
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change over cosmological times and hence can be assumed to be in equilibrium as a 
substructure. In any case, the morphology under discussion here will be ellipsoidal. 

Going from clusters to superclusters the following observation is worth noticing. 
Clusters are often elongated within their group [10] and aligned within their supercluster 
|41j . This filamentary network (-e.g. between the Virgo and A1367 cluster [IB]) 
could be a virialized substructure (including the clusters themselves) of the supercluster. 
Indeed, distinction between virialized and non-virialized regions of a supercluster exits 
|4*7j . In such a case these substructure would be best described by a prolate ellipsoid. 
If the angular momenta of the clusters are correlated [1H] we could have a model of 
rotating prolate/triaxial ellipsoid of clusters. 

Finally, even though superclusters are, for different reasons, not believed to be in 
equilibrium E3 ES] , the assumption of hydrostatic equilibrium still is often made in 
literature, either to test its grounds or to have a definite model [HH] • Even in cases where 
hydrostatic equilibrium is shown to be not valid, the results form the virial theorem seem 
to be still useful [HO] when applied to superclusters. This conclusion is corroborated in 



4. Major axis rotation: triaxial and oblate ellipsoids 

In this section we will consider the commonly known kind on ellipsoids configuration 
associated to cosmological and astrophysical systems, namely, oblate and triaxial 
configurations whose rotation is along the minor axis. Certainly a great portion of 
the number of cosmological structures have a non spherical symmetry, which can be 
related to internal streaming or rotation. 

Let us consider configurations with constant angular velocity along the minor axis. 
The kinetic energy tensor is then a contribution from internal motions and the rotating 
kinetic energy tensor, i.e. 

Tik — » T ik + lZ ik , (9) 

where lZik is defined for constant angular velocity as 

In the following analysis, we will neglect internal motion. The tensor virial equation 

for an homogeneous rotating configuration with angular velocity fl = Qe z axis is then 

8 c 
tt 2 vot (l ik - 5 iz l iz ) - \W ik \ + -rrp va , c 2ik = ~5 ik / Pd 3 r, (11) 

3 Jv 

where Wik is the gravitational potential energy tensor and lik is the moment of inertia 
tensor. From this expression one can derive two relevant equations. In the first place, 
we find the equation for the angular velocity, given by 

n2 \m x \-\m z \ , 8 /i zz \ 

^rot = + ^vp vac I 11, (12) 

which reduces to the standard Maclaurin formula (in the oblate case) for homogeneous 
configurations putting p vac = 0. On the other hand, by eliminating Q 2 from Eq.(fTTj). 
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we obtain a geometrical relation which restricts the geometry allowed for equilibrium 
configurations: 

T yy i\W xx \ - \W ZZ \) -X xx (\W yy \ - \W ZZ \) = (13) 

g 

~7rPvac-^zz (Ixx ^yy) • 

For a triaxial ellipsoid with constant density p and semi-axis a 1; 0,2 and 03 with 

ai > a 2 > a 3 , (14) 
Eq. f!13p can be recast into the following form 



2 = g 3 2 (^3(g2,93)-|C) (15) 

2 A 2 {q 2 , q 3 ) - A 1 (q 2 , q 3 ) + q 2 {A 3 {q 2 , q 3 ) - |C) ' 



where 



C = 2pvac/p, q i = a i /a 1 , (16) 

with q 2 > q 3 . The functions Ai are defined in dependence of q 2 and q 3 as (see |2HIE21ES]) 



A 2 (q 2 ,q 3 ) = 2g 2 <?3 



^3(^2,^3) = 2g 2 <?3 



'£(0, fc) - (1 - k 2 )F(9, k) - fk 2 sin 9 s 
(1 - k 2 )k 2 sin 3 9 t 
'^sin9 - E(9,k)\ 



[l -k 2 ) sin 3 9 /' 
with F(9, k) and E(9, k) are the incomplete elliptic integrals 

,2 



E(9,k) = { Jl-k 2 sin 2 
F(0,fc) = / % 1 zd0, 



?1 



1 - k 2 sin 2 V - 9a 

with 9 = arccos(g 3 ). Equation (fT3*|) is a complicated equation whose solution is 
a function g 2 = 92 (^3) representing the allowed ratios for configurations in virial 
equilibrium with major axis rotation. In figure Q we plot these solutions with and 
without cosmo logical constant. Note that the Maclaurin ellipsoid (q 2 = 1) is always 
a solution of (fT3|) . for any value of q 3 . Jacobi ellipsoids (triaxial case) are allowed 
for values of q 3 below the bifurcation point q 3 < q^ l{ . For A = 0, the bifurcation 
point is located at q 3 l{ = 0.583. On the other hand, with cosmological constant, we 
find the following relevant effects. First, the bifurcation point becomes a function of 
p vac . As long as we increase the value of £> the bifurcation point grows. This means 
that the separation between Jacobi and Maclaurin takes place at lower eccentricities. 
For instance, for ( = 0.6 we get q 3 l{ ~ 0.658. Another interesting effect is that the 
range of allowed q 3 also changes. Jacobi ellipsoids configurations are now restricted to 
q,mm < ^ < ^bif B e i ow g™ m ; the only allowed solution is again the Maclaurin ellipsoid. 
For ( = 0.6, the minimum value of q 3 is g™ 111 = q 2 ~ 0.24. Notice the value of ( is a time 
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Figure 2. Close plot for the second bifurcation point allowed for A ^ 0. 



dependent quantity, which depends on the density, p, of the astrophysical object we are 
investigating as well as on the value of p vac . Since we are interested in systems which 
have already reached an equilibrium configuration we consider ( as measured today. 
Hence, p vac = |pb g , where pb g is the background matter density. Taking into account 
the classical result from the spherical collapse model in an Einstein-de Sitter universe, 
then p « 200pb g . This gives ( ~ 0.023. Of course, this value may grow or decrease 
depending on the lower or higher density of the astrophysical system we are studying. 

In figure El we enlarged one plot for lower q% from figure 1 to explicitly show the 
second bifurcation point for £ = 0.1. For ( = 0.06, the second bifurcation point is located 
at q 3 ~ 0.12, while for smaller values of £, g™ m goes to zero as expected. For practical 
purposes, the behavior of the bifurcation point g!j can be approximately described as 
a function of £ by a fit. We find in good approximation 

qf(Q = 0.4826 (o.2082e 8 207C + l) , (19) 
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a 


b 


c 


qf 


„min 
% 





0.34622 


2.3035 


0.334876 


0.582 





0.006 


0.329658 


2.3639 


0.314573 


0.585 


0.026 


0.06 


0.332745 


2.3328 


0.318074 


0.589 


0.119 


0.2 


0.348911 


1.9253 


0.334721 


0.657 


0.335 



Table 1. Fits for the equilibrium condition (72(93) = ae bq3 — c for different values of £ 
for the triaxial-oblate solution. The range of validity of Jacobi's solution is also shown. 



n 2 n 2 

^"P ( = 25* P <T = 0.006 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

a 2 a 2 



2xp ( = 0.06 oil"'' (=0.2 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 3. Angular velocity of Maclaurin and Jacobi solutions for different ratios 
The black line corresponds to the Maclaurin solution. The red line is the Jacobi 
solution. 



while the equilibrium condition can be fitted to be approximately 

<? 2 (<? 3 ) = ae b " - c, (20) 

where the factors a, b and c (to be fitted) take different values as a function of the ratio 
(. In table d we summarize these values and our results on the bifurcation point of 
ellipsoids with major axis rotation. This bifurcation points separate triaxial solutions 
from oblate ones (a x = a 2 ). Finally, in terms of the quantities Aj defined at the 
beginning of this section, the angular velocity is written as 

Ss* = A 1 (q 2 , q 3 ) - qlA 3 (q 2 , q 3 ) - ^ U ~ |) , (21) 

In figure El we plot the angular velocity against the eccentricity 



<l - q\ (22) 

The Maclaurin and Jacobi solutions are shown. The explicit dependence of ( in the 
bifurcation point can also be seen. As in figure Q one can infer from this plot the 
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allowed ranges (in eccentricity) for the Jacobi's solution. One can also check the second 
effect of a non zero A on rotating configurations, which is to reduce the angular velocity 
with growing value of ( (see [Hp- 

To make contact with real properties of galaxy clusters we note that eccentricities 
of these objects have been extracted from simulations or observations at different red- 
shifts z j^SJ |Hj. In some scenarios the distribution of clusters around eccentricity 1 
is non-negligible [43] . If these objects are triaxial then, however, their eccentricity is 
limited to e(g™ in ) according to table 1. On the other hand, the change in q^ l{ will not 
be very significant. In view of the many effects of A on local astrophysical properties 
mentioned in the introduction, this is however apriori not obvious. 



5. Minor axis rotation: triaxial and prolate ellipsoids 

A note on convention is in order here. We define the prolate configuration by 

ai> a 2 = a 3 (23) 

to make a connection with the triaxial convention for semi-axis (fT4"|) i.e we would like to 
call in both cases ai the major axis. This is in contrast to the convention a\ = a 2 < a 3 
used in which, given (JHJ), is unsuitable to discuss a connection between prolate and 
triaxial ellipsoids as we will do below. Note that if we use the function Ai from reference 
j^Sj we have to interchange A\ «-> A 3 while keeping the eccentricity to be e 2 = 1 — q\. 
The prolate geometry offers a better scenario to test the effects of a positive cosmological 
constant. This can be seen by considering a non rotating homogeneous configuration 
|14j . By solving for the mean velocity from the kinetic energy we can write [Hj. 



<«>> - " 2|w 



2|l/wN| 

(24) 



2 \ 



2M 

This expression is valid for any geometry. The geometrical factor A defined as 
A = (16tt/3)(J/|VV) where VV = 2p~ 2 W and X = p' l l, is what enhances the effects 
of A: in the limit of high eccentricities, one can determine A Q u — > -^q^ 1 for oblate 
(q 3 <C 1) and A pro — > |g| [In (2q 3 )}~ 1 for prolate (q 3 ^> 1). The factor g| in found in the 
limit a 3 > a\ can enhance the effects of the density factor £, in contrast to the oblate 
case, where the ratio is only found to the first power. It is then expected that minor 
axis rotation which is not possible for an oblate spheroid might have some large effects. 

Minor axis rotation is a rare case, since it is expected that the flatness of a 
configuration is in part due to the centrifugal forces acting on a rotating configurations, 
as is the case of oblate systems. Nevertheless, such behavior has been observed in 
galaxies [23] and we can therefore expect to find a similar behaviors in galaxy clusters 
which are known to rotate JH]- Whatever the dynamical reason for such a rotation, it 
is interesting to investigate this issue from the equilibrium point of view. It has been 
stressed in the literature that a minor axis rotation is not possible for oblate geometries. 
However, it was not realized that the triaxial case is not an equilibrium solution unless 
we modify the virial equation, for instance, by the inclusion of A. 
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Figure 4. Bifurcation Prolate- Triaxial: the line q2 — 93 represents the prolate 
solution. 



The procedure of this section is similar to the previous section. We will explore the 
equilibrium condition and possible solutions in the case when the angular velocity lies 
along the major axis e x , in accordance with the convention chosen in [SHj for the triaxial 
ellipsoid (a\ > 02 > 0.3). The expression for the angular velocity is written in this case 
as 

o2 |ny - \w xx \ 8 / i xx \ 

fi rot = J o^Pvac I 1 - — J , (25) 

-L zz O \ -Lzz ' 

while the geometrical restriction (fTB^) becomes 

lyy (\m z \ - ) - X zz (\W yy \ - \W XX \) = (26) 

g 

-7t Pvac I xx (l ZZ -lyy), 

In terms of the functions defined in (|17|) . these expressions are written as 
3* = A 3 (q 2 , q 3 ) - qfA^ q 3 ) - |c (l - <h 2 

2 = q!(A 1 (q 2 ,q 3 )- fC) (2?) 

^1(92,93) + 93(^2(92,53) ~A 3 (q 2 , q 3 )) ~ |C 
Let us concentrate on the second line of Eq. (|2~7jl . Clearly, the prolate solution given 
in this case as q 2 = q 3 is trivially satisfied in with X zz = T yy and W zz = W yy , 
independent of the value of the ratio (. One also notices that in the standard case 
A = 0, the prolate solution is the only non trivial allowed figure in equilibrium with 
minor axis rotation. This implies that a triaxial configurations with minor axis rotation 
is not permissible as a configuration of equilibrium with zero cosmo logical constant. 

Switching on a non-zero cosmological constant changes the situation. Surprisingly 
for A 7^ we find a bifurcation point from which a triaxial solution can be determined. 
This is shown in Fig. EJ where we have plotted the solution q 2 (q 3 ) for different ratios (. 
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c 


qf 


„min 
% 


0.0001 


0.0024 


4.25 x 10" 5 


0.001 


0.00868 


4.25 x 10~ 4 


0.01 


0.0338 


0.00423 


0.1 


0.15 


0.0478 



Table 2. The two bifurcation points for the transition triaxial-prolate for different 
values of Q. 



si 2 a 2 

2xp f= 0.0001 \"P ( = o.ooi 




0.01 0.02 0.03 0.04^ 0.05 0.1 0.15 0.2 



Figure 5. Angular velocity of prolate and triaxial solutions for different ratios (. The 
black line corresponds to the prolate solution . The red line is the triaxial solution 
allowed only with cosmological constant. 

In this figure, the straight line represents the prolate solution q% = q2- The bifurcation 
point arises for small values of q 3 , which represents large values in the eccentricity 
e = \A — <?i- F° r the largest value of £, £ = 0.1, we obtain the largest value for the 
bifurcation point, located at ^3 m 0.15. Hence the triaxial solution arises for e > 0.98. 
Note that the triaxial solution has also a lower limit of validity for q 3 . It occurs when 
q2 approaches 1. For instance, for ( = 0.1, the triaxial solution is valid in the interval 
0.047 < q 3 < 0.15 (0.988 < e < 0.998). As long as ( decreases, the range of validity of 
the triaxial solution is reduced. The complete set of values of our solutions are given in 
table El 

In figure El we plotted the angular velocity against q 3 . The bifurcation point, q 3 lf , 
is where the prolate curve intersects the corresponding triaxial case, g™ 111 is represented 
by the left vertical line beyond which we have no solution. It can be nicely seen that 
with decreasing ( the range of the triaxial solution becomes smaller until it vanishes 
for C = 0. The angular velocity is also not a smooth function as it was the case for 
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the triaxial/oblate case. Indeed, it oscillates (allowing for the moment the physically 
excluded range of Q < 0) when approaching g™ in . As can be seen from figure 03 the 
restriction of positivity on Q narrows the allowed range of ^3 further. Given the fact 
that without A triaxial configurations with minor axis rotation are not in equilibrium, 
this narrow range is rather a virtue than a disadvantage if we want to put the theory 
to test. Had we an accurate knowledge of a triaxial ellipsoidal cluster rotating the 
minor axis, the theory outlined above could be easily verified or falsified, of course 
bearing in mind the simplified assumptions of constant density/no internal motion (as 
in [2E]) and the fact that only projected information of the axis is available to us. From 
the observational side, its encouraging to know that dark matter halos around galaxy 
clusters are described by triaxial configurations which tend to be more prolate |44j . 

6. Conclusions 

We have investigated effects of a positive cosmological constant on astrophysical objects 
with ellipsoidal geometries. We have considered both commonly known homogeneous 
ellipsoid configurations associated to cosmological and astrophysical systems, namely, 
oblate and triaxial configurations whose rotation is along the minor axis, as well as 
configurations with constant angular velocity along the minor axis. 
We made use of the Newtonian limit of the Einstein field equations and have solved the 
tensor virial equation for the angular velocity. 

Our results generalize the ones obtained in [2E| which applies for the case A = 0. This 
is especially true for the bifurcation points of rotating ellipsoids as solutions of tensor 
virial equations for the oblate geometry. Adding to that we have also investigated the 
prolate/triaxial case with minor axis rotation, which has not been examined before as 
this case is only possible with A / 0. 

We have found that the bifurcation point, which in the oblate case distinguishes the 
Maclaurin ellipsoid (oblate) from the Jacobi ellipsoid (triaxial), is sensitive to the 
cosmological constant. A new effect is given by the appearance of g™ m which restricts the 
validity of the solutions from below. This point exists only for a positive cosmological 
constant. We also show that the cosmological constant allows triaxial configurations 
of equilibrium rotating along the minor axis as solutions of the virial equations. The 
significance of the result lies in the fact that minor axis rotation is indeed found in 
nature. Being impossible for the oblate case, it is permissible for prolate geometries, 
with A zero and positive. However, for the triaxial case an equilibrium solution for the 
minor axis rotation is found only for non-zero positive A. Therefore the latter case is 
intimately, albeit indirectly, related to the existence of a cosmological constant. 

The ellipsoidal geometry which we have studied in this paper is well suited to 
model many geometries, such as the disk. Therefore, the results of this article should 
in principle be applicable also to superclusters which are usually found to be very flat 
objects [19]. 
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